ONE NON-RELATIVISTIC PARTICLE COUPLED TO A 

PHOTON FIELD 
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Abstract. We investigate the ground state energy of an electron cou- 
pled to a photon field. First, we regard the self-energy of a "free" elec- 
. tron, which we describe by the Pauli-Fierz Hamiltonian. We show that, 

in the case of small values of the coupling constant a, the leading order 
Mh I term is represented by 2a7r _1 (A — m[l + A]). 

f"^ ■ Next we put the electron in the field of an arbitrary external potential 

£SJ ' V, such that the corresponding Schrodinger operator p 2 + V has at least 

one eigenvalue, and show that by coupling to the radiation field the 
CO ■ binding energy increases, at least for small enough values of the coupling 

constant a. Moreover, we provide concrete numbers for a, the ultraviolet 
cut-off A, and the radiative correction for which our procedure works. 
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In recent times great effort was made to start a program with the goal 
Q^. to put Quantum Electrodynamics (QED) on a firm mathematical footing. 

Bach, Frohlich, and Sigal were the first who proved (in |BFSUlEFS51lBFS3] ). 
at least for small values of various parameters, that atoms and molecules 
have a ground state in the presence of a quantized radiation field. In the 
one particle case Griesemer, Lieb, and Loss in |GLL| succeeded in removing 
the restrictions on the parameters. We refer to [G"L~L~] for an extensive list 
of references concerning this subject, 
jjj ■ In the present paper we are especially interested in the behavior of one 

electron coupled to a photon field. We consider the "free" case as well as 
the presence of an arbitrary external potential V with at least one negative 
eigenvalue in the electron case. We are interested in the ground state energy 
of such an electron when the coupling constant a is small and the ultraviolet 
cutoff parameter A is fixed. 

Without radiation field the self-energy of an electron is simply given by the 
bottom of the spectrum of —A, which is equal to 0. In the case of coupling 
to a photon field the situation changes dramatically and the self-energy is a 
complicated function depending on a and A. The fact that evaluating this 
quantity is a highly non-trivial problem was pointed out in |LL| by Lieb and 
Loss, who, in contrary to us, considered the case of fixed a and large values 
of the cutoff parameter A. 
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The paper is organized as follows. In the present section we state our 
main results on the self-energy and on the binding energy of an electron, 
which will be proved in Sections 3 and 4. In Section 2 we comment on our 
notation and in Section 5 we evaluate the radiative correction in the case of a 
nucleus. We provide concrete numbers for a and A for which our procedure 
works. 

1.1. Self-energy of an electron. The self-energy of an electron is de- 
scribed by the Pauli-Fierz operator 

T = (p + \faA(x)) 2 + \fa,G~B{x) + Hf. (1.1) 

We fix units such that H = c = 1 and the electron mass m = \. The 
electron charge is then given by e = ^fa, with a ~ 1/137 the fine structure 
constant. In the present paper a plays the role of a small, dimensionless 
number which measures the coupling to the radiation field. Our results 
hold for sufficiently small values of a. a is the vector of Pauli matrices 
(°i> 02) 03)- Recall, the (Tj are hermitian 2x2 matrices and fulfill the anti- 
commutation relations OiOj + OjOi = 2I C 28ij. The operator p = —iV is the 
electron momentum while A is the magnetic vector potential. The magnetic 
field is B = curl A. 

The underlying Hilbert space is 

H = £ 2 (K 3 ;C 2 ) ®-F 6 (£ 2 (IR 3 ;C 2 )) (1.2) 

where Th is the bosonic Fock space for the photon field. 
The vector potential is 

A— 1,2 

and the corresponding magnetic field reads 

b( x ) = E t L lW {k A Z£x{k)) [ ax{kykx " a *^ e ~ lkx ] dk > ( L4 ) 

where the operators a\ and a* x satisfy the usual commutation relations 

[a u {k),a* x (q)] = 5(k - q)5 x ,v, [a\(k), a v (q)] = 0, etc. (1.5) 

The vectors £\(k) £ M 3 are any two possible orthonormal polarization vec- 
tors perpendicular to k which we assume to be either odd or even, i.e. 

e x (-k) = ±e x (k). (1.6) 

The function xd&l) describes the ultraviolet cutoff on the wavenumbers k. 
We choose x to be the Heaviside function 0(A— \k\/lc), where lc = h/(mc) 
is the Compton wavelength. Therefore A = 1/4 corresponds to the photon 
energy mc 2 , which sometimes is considered as a natural upper bound for 
the cut-off parameter, since it is the maximum value that guarantees that 
no pair-production takes place. 



ONE NON-RELATIVISTIC PARTICLE 



3 



The photon field energy Hf is given by 

H f =J2 f \k\a* x (k)a x {k)dk. 

We estimate the leading order in a of the ground state 

E Q = inf spec T. 

Recall that a{0\A 2 \0) = avr^A 2 , where |0) is the vacuum in J-j,, describes 
the vacuum fluctuation of the field A. Thus, this term is somehow ab initio 
included in the self-energy. Therefore it may turn out as a slight surprise 
that a7r _1 A 2 is not the leading order term of £ a , which relies on the fact that 
the magnetic field yfaaB reproduces the same number but with a negative 
sign. We show that the leading order term is given by 2avr- 1 [A-ln(l + A)], 
which is the content of our first theorem. 

THEOREM 1. Let A be fixed but arbitrary. Then 

|S Q - 2a7r" 1 [A - ln(l + A)]| < da 2 {A 2 + A 3 ) + C 2 a 3 (A 3 + A 4 ) ln(l + A) 

(1.7) 

for some constants C\i > independent of a and A. 

It turns out in the proof that one photon is enough to recover the leading 
order term in a. 

1.2. Increase of the binding energy. We consider an electron in an arbi- 
trary external potential, described by a real valued function V ', independent 
of a, which satisfies (i) that the negative part of V, V_, is dominated by the 
kinetic energy p 2 = —A, 

y_ < -eA + C £ , (1.8) 

for any e > and for some positive constant C e , such that the self-adjointness 
and boundedness from below of p 2 + V is guaranteed. Additionally, since 
we want the essential spectrum to start at 0, we require (ii) that V- tends 
to at infinity. Furthermore, we assume (iii) that the operator p 2 + V has 
at least one negative energy bound state. 

The corresponding operator with a quantized radiation field reads 

H a = T + V, (1.9) 

for which |Hij guarantees self-adjointness. Recently, Griesemer, Lieb, and 
Loss jGLLj proved for all values of a, A that H a has an eigenstate, as well 
as that the binding energy S a — E Q , with 

E Q = inf spec H a , (1-10) 

cannot decrease, i.e. 

E^-S^-eo, (1.11) 
where — eo = inf spec \p 2 + V]. We show that at least for small values of the 
coupling constant a the binding energy increases, namely (|1.11|) holds with 
strict inequality. 
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THEOREM 2. Let V satisfy the conditions (i)-(iii), and <f> be the ground 
state of p 2 + V with corresponding energy — eo- Then 

B a -T, a < -eo - aS(y, A) + 0(a 2 ) (1.12) 

for some positive number £ (V, A) depending on V and the cut-off A. More- 
over, there exists a number p > 0, such that for all a £ (0, p] 

E a - S a < -eo. (1.13) 

REMARK 1. The equation (fTT^ holds for all real V, such that p 2 + V 
is self-adjoint and has a ground state. But only for V satisfying (i) — (iii) 
we know that H a has a ground state. 

We actually check that S(V, A) is bounded from above by 

£(V,A)<\\p<t>\\ 2 ^Hl + A]. (1.14) 

In the spinless case (B = 0) it was previously proven ( |HWj ) that cou- 
pling to a photon field can produce a ground state even if the corresponding 
Schrodinger operator has only continuous spectrum. This phenomenon of 
enhanced binding was earlier shown in |HSj in the case of dipole approx- 
imation (A = ^4(0)) in the limit of large values of the coupling constant 
a. 



2. Notation 

Throughout the paper we use the notation 

A(x) = D{x) + D*(x), B(x) = E{x) + E*{x) (2.1) 

for the vector potential, respectively the magnetic field. 

The operators D* and E* create a photon with wave functions G(k)e~ lk ' x 
and H(k)e~ lk x , respectively where G(k) = (G 1 (k) , G 2 (k)) and H(k) = 
(H 1 (k) , H 2 (k)) are vectors of one-photon states, given by 

and 

HX{k) = i^ kA£x{k) - (2 - 3) 

Applying Fourier transform, with respect to k, on G(k)e~ lk ' x , i.e. 

F- x [G{k)e- ih - x ](£) = [ e ik <t- x) G(k)dk = D*\0)(£-x), (2.4) 

we see that D* creates a photon depending on £ — x, where £ denotes the 
photon variable in configuration space. 

If the reader is not used in thinking of a "position" of a photon, he may 
consider it as a simple unitary transform. Because, now we introduce the 
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variable r] = £ — x and the object we are going to work with is the Fourier 
transform with respect to 77, namely 

F[D*\0)( V )](k)=G(k). (2.5) 

The same result holds for E* if G(k) is replaced by H{k). 

Notice that in the so chosen notation the electron momentum p = —iV x 
of a one photon state D*f(x,rj) = f(x) (g) D*\0)(rj), with 77 = £ — x, can be 
written as 

3 

pD*f(x, V ) = (P*/0»0) ® A*|0)(r?) - /(x) 8) {p v D*\0)( V )) 

i=l 

= \p x - Pv ]D*f(x, V ). (2.6) 

It turns out to be convenient to denote a general vector £ Ti as a 
sequence 

* = {^V'l,^,-}, (2-7) 

where 

V> n = i^n(x,Ci -x, ...,£ n -x), (2.8) 

with relative variables £j — x = r/j. For simplicity, we suppress the polariza- 
tion of the photons and the spin of the electron. 

3. Proof of Theorem [T] 

3.1. Upper bound. We choose the sequence of trial wave functions 

*n = {fn(x) T, -v^K + H f ]-\a ])E*f n (x, 77), 0, 0, . . .}, (3.1) 

with 77 = £ — x, f n 6 £ 2 (M 3 ), 1 1 / n || = 1 and | denoting the spin- up vector 
in C 2 . We assume for the electron part f n {x) that ||p/ n || — ► as n —> 00, 
which means that the electron spreads out to infinity. This is reasonable, 
for in the case without coupling this sequence recovers the bottom of the 
spectrum of —A. Recall, 

E*f n (x, V ) = f n (x)®[E*\0)]( V ) (3.2) 

and, due to 77 = £ — x, 

(b 2 + H f ]E*f n ) (x, 77) = (\p x - Pv } 2 + \p v \)E*f n (x, 77). (3.3) 

Notice, in the one photon case Hf = \pr,\. Observe, 

[p,D}= V f G x (k)-[p,e ikx ]a x (k)dk = 

= Y [ (G x (k)-k)e ik - x a x (k)dk = (3.4) 

X=l,2 

by our choice of e\(k). Analogously, 

\p,D*] = \p,E] = \p,E*}=0. (3.5) 
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By means of Schwarz' inequality it is easy to see that by our assumption 
||p/n|H0 

lim (^ n ,pA^ n ) = 2 lim R(¥ n ,p£>tt n ) = 0. (3.6) 

n— >oo n— >oo 

Similarly, if we denote the 1-photon part of \& as 

= -^\pl + H f \-\a ])E*f n , (3.7) 

then 

= \p x + Pv ] 2 W) = m,p 2 r,W) + Errov(p x f n ), (3.8) 
with limn^oo Error(pa;/ n ) = 0. Furthermore, we compute 

($S [p\ + Hf]tf) + 2^{(a })E*f n ,^) 

= -a((a ])E*f n , [p 2 + Hf]~ 1 (a ])E*f n ) 

= -a(0\(a ])E[p 2 + H f ]-\a ])E*\0) (3.9) 

In momentum representation, recall J-[E*\Q)]{k) = H(k), the last term in 
(|3.9j) equals 

(0\(*VE[P 2 + H f ]-\*VE*\0) = P^^^ dk 

= I WTW\ dk = imp2 + H ^ lE * l0) - (3 - 10) 

The equality from first to second line holds, because of the anti-commutation 
relations of a and the fact that H is purely imaginary. 
For simplicity, we will use throughout the paper 

\H(k)\ 2 = ]T \H x (k)\ 2 , \G(k)\ 2 = £ |G A (A0| 2 . (3.11) 

A=l,2 A=l,2 

We evaluate 

<o|^ + ^,-^|o>=E 2 /J^M* 

= sj/K* - 2 ' A - ta(1 + A)11 ' (3 ' 12> 

Since 

[Z>, £)*] = DD* - D*D = vr^A 2 , (3.13) 

we know 

A 2 = (D + D*) 2 = 7r _1 A 2 + 2D*D + D*D* + DC. (3.14) 
With (^ n ,DD^ n ) = and lim„_ +00 (?/;'f , ^") < 2avr _1 ln(l + A) we arrive at 

lim (¥„,T¥ n ) < avr" 1 A 2 -a(0|£;[p 2 + iJ / ']" 1 J B*|0) + 2vr" 2 a 2 A 2 ln(l + A) 

n— >oo 

+ 2a 2 || J D[p 2 + i7 / ]~ 1 (cj |)£*|0)|| 2 . (3.15) 



ONE NON-RELATIVISTIC PARTICLE 



The last term in the r.h.s. vanishes which can be seen by explicit calculations 
using the relation 

E 44 = *„■ - (3-16) 

A=l,2 1 1 

Namely, if e J ' n denotes the totally antisymmetric epsilon-tensor, we obtain 
f G>>{k)Hj{k) 3 . f x (\k\)s\(k)[eJ^s[(k)k n ] 

» rx(\k\)[Si,i-m^ ln k n 
= y i ,,,, '* dk = o. 3.i7 

Therefore, since (^ n , = 1 + const. a, we obtain the upper bound 
S Q < lim(¥ n ,r* n )/(tf ni tf n )< lim(tf n ,Ttf n )< 

n^oo n— >oo 

< 2avr~ 1 [A - ln(l + A)] + 2tt" 2 a 2 A 2 ln(l + A). (3.18) 

3.2. Lower bound. We have learned in the previous section that an ap- 
proximate ground state ^f, \\^f\\ = 1, satisfies 

TV) < 2avr~ 1 A + 2yr~ 2 a 2 A 2 ln(l + A). (3.19) 

Since 

T= [a- {p + V^A(x))} 2 +H f > H f (3.20) 

we immediately get 

(*, H f V) < 2air~ 1 A + 27r~ 2 a 2 A 2 ln(l + A). (3.21) 
Therefore, by means of Schwarz' inequality we obtain 

2y/a(V,pAV) = 4y/a$t(y,pDV) < a\\p^\\ 2 + 4a~ 1 a||L>1'|| 2 (3.22) 
y/a(^,aB^) = 2^(^,aE^) <cA 2 a\\^\\ 2 + c~ 1 (l/A 2 )\\E^\\ 2 , 

(3.23) 

for any a, c > 0. Since, by A 2 > 0, 

llp^ll 2 < {V,Tty)-{V,H f y)-4:y/a&(y,pDy)-2y/a(it>,aEy), (3.24) 

we obtain by ()3.22|) . (|3.23|) . the operator inequalities (e.g. |GLLl Lemma A. 
4]) 

D*D < -Ail/, E*E < 3^A 3 iJ/, (3.25) 

and (l3~2TT> 

\\p^\\ 2 < C x a(A + A 2 ) + C 2 a 2 (A 2 + A 3 ) ln(l + A) (3.26) 

for some suitable constants C± t 2 > 0. Equations ()3.21|) and (|3.26j) will be 
decisive to control our error estimates. Recall, 

A 2 = tt~ 1 K 2 + 2D*D + 2^kDD. (3.27) 
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Thus, since D*D > 0, 

oo 

Ttf) > avr-^ 2 !!*!! 2 + f [^o, H + X^n, V>n+i, ^n+2], (3.28) 

n=0 

where 

£obh,H = {tpi,Atpi) + 2V^maE* + 2pD*]iPo,tp 1 ) (3.29) 

and 

+ 25ft (([yfeaE* + 2^pD*]^ n+1 + aD*D*^ , Vn+2) , (3.30) 

with 

A=p 2 + H f . (3.31) 
Recall, we can write our approximate ground state \P as 

* = {ipo(x),il)i(x,rn), ....,ip n (x,rix> •••>%)> •••}> ( 3 -32) 
with 77j = — x. For convenience we will work in momentum representation 

^[ip^X,^!, ...,7] n )](l,ki, ...,kn) = 1p n (l,ki,.. . ,k n ), (3.33) 

with k = (k\,...,k n ) and I denotes the momentum representation of the 
electron variable. 

We consider the term £ [i(} n , ip n +l, ifin+z]' As a straightforward conse- 
quence of Schwarz' inequality we derive 

£[ip n ,ip n+1 ,ip n+ 2} > 

>- s/a[A" l/2 aE* + A~ 1/2 pD*]t/j n+1 + aA- l/2 D*D*^ n ' \ (3.34) 
Similarly, 

£c#o,^i] > -a\\(A- l/2 aE* + 2A- 1/2 pD*)^ \f (3.35) 

The fact that H(k) is purely imaginary and A commutes with the reflection 
I — ► — / imply 

5(E*ip ,A- 1 E*Tp )=0. (3.36) 
Together with the anti-commutation relations for a we infer 

(crE^A^aE^o) = (^EA^E*^). (3.37) 

Now we are going to evaluate the r.h.s. of (|3.34|) . These evaluations then 
are also applied to (|3.35j) . First, we consider the diagonal terms. The most 
important one, that reproduces the leading order in a, is 

-a((rE*tlj n+1 ,A~ 1 (rE*i; n+1 ). (3.38) 
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Explicitly, in Fourier representation, 

, n+2 

F[E*ip n+1 ](l, h,..., k n+2 ) = , t , y~) H(ki)ip n+1 (l, fci, , k n+2 ), 

(3.39) 

where tyi indicates that the i-th variable is omitted in tp n +i- Due to per- 
mutational symmetry the expression (|3.38f) consists of two different terms 



-a(aE*il; n+1 ,A- l aE*il> n+1 ) = -a(l n+1 + II n+1 ). (3.40) 

The diagonal term I n +\ where in the right as well as left hand side of ()3.38j) 
a photon is created with index i and the mixed term II n +i- By similar 
arguments as for (|3.35|) . namely if we take advantage of the fact that H is 
purely imaginary and use the anti-commutation relations of a we obtain 

W < / ^(W 2 )| 2 y rl (/, fel ,- + ^ +1 )P 

i [z-ES ^] 2 + E?=i N 

We set Q = [I- ES. 1 + 1 ^ +2 | 2 + E^ N and 6 = 2[l - Zti *»] • 
and expand 

1 1 1,1 1,1,1 , , , 

= — + —b— + —b b—. (3.42) 

Q-6 Q Q Q Q Q-b Q V ; 

Plugging this expansion into (j3.41j) . the second term, when integrating over 
k n +2, obviously vanishes. Since Q > |fc„+2| 2 + |&n+2| and Q — b > \k n+2 \ we 
estimate 

/„ + i<||Vwi|| 2 f Mr^jdk 



\k\ 2 + \k\ 

+ 4 f \H(k n+2 )\ 2 \k n+2 \ 2 \l - Zti kj\ 2 \jJn+i(l, ku • • • , K +1 )\ 2 
J [\k n+2 \ 2 + \k n+2 \] 2 \k n+2 \ 

< H^n+ill^Ol^-^lO) + const.A||pVn+i|| 2 . (3.43) 
For the mixed term we have 

JT <K, + llV [J^WMLv 

IW» < 3(n + l) £j _ g? __ gf _ x 

x \ip n+1 (l,k 2 , . . . , fc n+2 )||V'n+l(^ fei, ■ ■ ■ , k n+1 )\dldki . . . dk n+2 . (3.44) 
By means of the one-photon density 

Pi> n+1 {k) = (n + 1) J \ip n+1 (l,k,k 2 , . . . ,k n+1 )\ 2 dldk 2 . . .dk n+1 (3.45) 
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we infer after applying Schwarz' inequality to (j3.44|) 



Hn+l <Hn + l)^2 [ 

\ 1 o " 



1^(^)11^(^+2)1 



A=l,2 • 



< 



< 



\ki\ + |fc n+2 | 

x \J P^n+i \J Pi>u+i {k n +2)dkidk n+2 
6(n + l) f x(\h\)x(\k n+ 2\) A"- TTT 

~WrJ \h\ + \k n+2 \ y/M<>**M* 

x yj | k n+2 1 Pijj n+1 {k n+2 ) dki dk n+2 

mm^ dkldkn+2 



^2 / \ k \Pi>n+l( k ) dk 



1/2 



(N + \kn+2\r 

< const.A 2 (V„+i,^/Vn+i). (3.46) 



Therefore, we summarize 



a{aE*^ n+1 ,A- l aE*i> n+ i) > -a{0\EA~ l E* \v, \\ Wn+1 \ 

- const.aAf ||p^„+l|| + A(ip n+ i, H f ip n+ i)) . (3.47) 



All other terms in (|3,34|) are of order 0(a 2 ) or even of higher order. 
For the second diagonal term we obtain 

-a 2 (ip n ,DDA- 1 D*D*ip n ) > -const.a 2 A(A||^ n || 2 + ||^ n || 2 + (^n,^/^ 

(3.48) 

The proof of (|3.48|) . whose strategy is similar to the one for (|3.47[) . will be 
postponed to Lemma the appendix. 
The third diagonal term of l|3.34j) reads 



-ax 



A=l,2 

G x (k n+2 )-(i-Y,tiki)] 2 \ipn+i(i,h, 



n+2 



-dldk\ . . . dk- 



n+2 



+ (n + 1) 









G A (fc n+2 )-(/-Er=i 2 ^ 






i - zr=i 2 h 


2 +£r=i 2 N 



X 1p n+ i(l, h,..., kn+l)lpn+l(l, k 2, ■ ■ ■ , k n+2 )dldki . . . dk n+2 



> — const. aA 



n+1 1 



+ (W i +i,#/Vwi))- (3.49) 



ONE NON-RELATIVISTIC PARTICLE 

For the second term in the r.h.s. we used first 

< 1, 
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2 


i - EIL + i 2 h 





(3.50) 



and afterwards applied Schwarz' inequality to bound it from below by 
> -(n + 1) J \° X ( k i)\ \ G \ k n+2)\ y /p^ n+1 (k 1 )p^ n+1 (k n +2)dk 1 dk. 

x(\h\)x(\k n+2 \) 



n+2 



A=l,2 ' 



> -const. / \k\p^ n+1 (k)dk 



n l/2 



"4 /t n+2 



-dkidk n+ 2 



, (3.51) 



which yields the last term in the last line of (|3.49[) . 
Next we consider the off-diagonal terms. Notice, 



aU(pD*ip n+1 ,A l aE*^j n+1 



A=l,2 



G x (k 


n+2) ■ {I ~ J2?=l fa) 


H x (k n+2 y 




i - £IL + i 2 h 


2 + Y2=x\h\ 



x -(V>n+i,cnAn+i)c 2 (^i> . . • , k n+ i)dldki . ..dk n+2 
+ (n + l)5R 





H\k n+2 y 




i - EIL + i 2 h 


2 + ^\ki\ 



x ■(il> n+ i(l,k2, ■ ■ ■ ,k n+2 ),cnpn+i(l,ki, . . . , k n+ i)) C 2dldki . ..dk n+2 . (3.52) 



The first term in the r.h.s. vanishes, because the integral is purely imaginary. 
For the second term we use J^_ b < ^6 -1 / 2 and Schwarz' inequality to bound 
it from above by 



(13321 < a 



|G A (fci)||tf A (fc n+2 )| 



? 2 J (\k 1 \ + \k n+2 \y/ 2 \k 1 mk n+2 \^ 



x yj PVn+i (fa ) I fc i I \J Pi^n+i {kn+2)\ k n+2 \dk!dk n+2 

< const.aA 3 / 2 (^ n+1 ,ff / ^ n+ i). (3.53) 

The second off-diagonal term, 2a 3 / 2 5ft(pD*VVt+i, A~ 1 D*D*tp n ), can sim- 
ply be bounded from above by 



a{^ n+1 ,pDA- l pD*^ n+1 ) + a 2 {^ n ,DDA- l D*D*^ 
on which we now apply Ij3.48j) and (j3.49j) . 



(3.54) 
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Finally, by means of Lemma El in the appendix we estimate the last off- 
diagonal term by 

a 3 ^(aE*^ n+1 ,A- 1 D*D*^ n ) < const.aA(a||Vn+i || 2 + (V>„, tf/Vn)) • 

(3.55) 

Collecting above estimates and summing over all n we arrive at 

(¥,Ttf) > 2cm~ 1 (l -const.a(l + A))A 2 ||^|| 2 - a||*|| 2 (0|E^l -1 ^*|0) 

- const. aA(l + A)(tf , HfV) - const. aA||p^|| 2 . (3.56) 
By our a priori estimates (|3.21f) and (|3.26|) we prove the theorem. 

4. Proof of Theorem E] 

It suffices to use a cleverly chosen trial wave function. We assume that 
<t>(x) G C 2 (M. 3 ) is the ground state of p 2 + V, i.e. 

(p 2 + V> = - eo 0. (4.1) 

Recall, 

3 

E*</>(x,r,) = 4>(x) ® [E*(x)|0)](77), D*p(f> = ^p^(x) ® [A*|0)](r?) 

i=l 

are 1-photon functions depending on x and the relative coordinates ^ — x = 
■q. Therefore, in configuration space, where if)(x,rj) denotes one of these 
functions, we have (cf. the previous section) 

[(p 2 + H f )iP](x,7]) = (\p x - Pv } 2 + |p„|)#c,77). (4.2) 

For sake of convenience we define the operator 

Av = vl +P* + H f + V + e = (p 2 + V + e ) ® I + I ® (pj + fl/), (4.3) 

which acts on £ 2 ((R 3 ; C 2 ) (8) (M 3 ; C 2 )) and is obviously positive and invert- 
ible. 

Now, we choose our trial wave function \& G 7i as 

$ = {0( x ), -2^A- 1 L»*p0(x, ??) - V^A^aE^ix, rj), 0, 0, ...}, (4.4) 

where <fi = 4> ] and rj = ^ — x. We assume ||0|| = 1 and for simplicity we 
denote the 1-photon part of *S as ip±. Notice, 
oE*4>(x,r]) = 4>{x) ® (<r T)_-E*|0)(?7). 
First, observe that Ay(f) = yields 

A v aE*4> = (p 2 + H f )aE*4>. (4.5) 
Therefore, since Ay and p 2 + iJj commute by definition, we infer 

Ay X aE*4> = (p 2 + HfY^oE*^ (4.6) 

and 

§,GEAy l aE*$) = U\\ 2 (0\E[p 2 v +H f r 1 E*\0). (4.7) 
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Thus, we evaluate 

(*,H a *) = avr _1 A 2 ||^|| 2 - e ||^|| 2 + (ipi,A v ipi) 

+ 2^U {{ctE* + 2D*p)4>, Vi) + 2a(^D*D^x) + 2(ik,p x ■ p^x). (4.8) 

Notice, the cross term 

aEA- l D* P 4>) = aE[p 2 + H f }- 1 D*p0) (4.9) 

vanishes as in the previous section, because it is purely imaginary. Therefore, 
by our choice of tpi we get 

(V?i, Avfa) + 2^ ((<tE* + 2D*p)0, ^i) 

= -4a(D*p x ^Ay 1 D*p x ^) - aU\\ 2 (0\E[p 2 v + Hf^E^O). (4.10) 
Moreover, we have 

(.i>l,Px -Prj^i) = a(D*p x cj),Ay 1 p x ■ p r] A v 1 D*p x (j)) 

- a(aE*0, \p\ + Hf}-^ ■ p n \p\ + Hf\~ l aE*4>) 

+ 2aM(D*p x $, Ay l p x ■ p n Ay X oE*4>). (4.11) 

The first term in the r.h.s. vanishes by integrating over the 77-variable (the 
best way to see it is using the representation in momentum space), due to 
(|1.6[) and the fact that the operator Ay commutes with reflection rj — > — 77, 
respectively k — ► —A;. The second term vanishes when integrating over the 
^-variable (notice, (<j>,p x <p) = 0), and the third term vanishes, because it is 
again purely imaginary. 

Furthermore, since DAy X aE*4> = (see (JSHZJ)), we realize, after straight- 
forward calculation, that 

a{if) 1 ,D*Dip 1 ) < const. a 2 ln(l + A) 2 . (4.12) 

Using above estimates and the fact that 

air^A 2 - a(0\E[p 2 + Hf^E^O) = 2a7r" 1 [A - ln(l + A)] (4.13) 

we infer 

H a *)/(*, *) < -eo + 2avr- 1 [A - ln(l + A)] 

- 4a(D*p x <P, A v 1 D*p x <P) + 0(a 2 ) ln(l + A)(A 2 + ln(l + A)) (4.14) 

which proves the first statement of the theorem with 

£(V, A) = ±a{D*p x <t>, Ay 1 D* Px( j ) ). (4.15) 

The second statement follows by the observation that (D*p x (f>, Ay 1 D*p x cf>) 
is strictly larger than 0, which is a consequence of the fact that D*p x cj) is a 
not identically vanishing function £ £ 2 (R 3 ® (M 3 ,C 2 )) and Ay an invertible 
operator. 
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5. Computation of concrete numbers 



5.1. Error for the self-energy. We are going to calculate the error Err 
of the self-energy, 



|S a -2avr~ 1 [A-ln(l + A)]| < £rr(a 2 



(5.1) 



First, it is important to estimate the kinetic energy term \\p^f\\ for an ap- 
proximate ground state with quite good constants. Since we will calculate 
anyway with values A < 1 we can assume ('J/, T^S) < 2an~ 1 A. 
By means of (|3.22j) . (|3.23j) and then applying (|H.25j) we obtain 

2aTT' 1 A > (¥,Ttf) > (1 - a)||^|| 2 + 

- 4a~ 1 a||L>1'|| 2 - caA 3 - l/(cA 3 )||£7tf || 2 



> (1 — o)||p*H 2 — cotA 3 + 



1 - a _1 87r _1 aA - c" 1 — 
3vr 



(*,H f V), (5.2) 



where 1 > a > 0. We require the last term [..] to be > 0. For simplicity, we 
choose c = -, then our first condition on a reads 



air 
4A 



and additionally 



WpW < 



2 ^ 2aA(l + A 2 ) 



7r(l - a) 



(5.3) 



(5.4) 



The main contribution to Err stems from the the terms (|3.47j) and (|3.49|) . 
In fact the the third diagonal term (|3.48|) is negligible compared to (|3.47|) 
and (JH30J). 

Evaluating the corresponding integrals yields 



MJ) < 2a7r~ 1 A||j9*|| 2 + 2a7r 1 A(*, #/<£), 



and 



(I3~i3]l + lEORjl < — aA||ptf|| 2 + — aA 2 (tf,fl/tf). 



(5.5) 



(5.6) 



In order to incorporate the off-diagonal terms we, for simplicity, double the 
so gained value for Err and derive 



Err < 2 

< a 

where we used (|3.21[) and (|3.26j) 



14 7 
— aA||^|| 2 + — aA 2 (*,H f *) 
6ir 2ir 



56A 2 (1 + A 2 ) _ 2a2 

+ 14-7T A 

3tt 2 (1 -a) + i47rA 



(5.7) 
(5.8) 
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5.2. Radiative correction. We consider an electron in the field of a nu- 
cleus with charge Z, i.e. 

\x\ 

where (5 is the "real" fine structure constant = 1/137. The ground state 
energy of the corresponding Schrodinger operator p 2 + V is known to be 

inf spec \p 2 + V] = -eo = ~\{(3Z) 2 . 
The corresponding radiative correction, obtained in Q4.15|) . is given by 

S(V,A) = -aA{4) :P DA- l pD*(t)), (5.9) 
where 4> = denotes the ground state of p 2 + V. We know 

V<f>(\x\) = d r 0(r)e r (0, 99) = el /2 <t>{r)e r {0, <p) (5.10) 

when using polar coordinates. Denote with ^ the eigenstate of p 2 + V with 
corresponding eigenvalue — e,. Then, by means of (|5.1U|) . we obtain, by 
straightforward computations, 

(frpDAfpD*® > 4ire V |q| 2 / P 2 dp = e F(A), (5.11) 

jXL Jo e o - ei + + p 

where 



&(r, 0, ¥?)^(r) cos{9)r A drdVL, (5.12) 

with = sin8d8dip, and the sum runs over all hydrogen eigenstates. No- 
tice, due to textbooks ^ |q| 2 ~ 2/15 = which in physicists' words is 
expressed by "80 percent of the average excitation energy (eo — e n )AV is 
achieved by the continuous spectrum". Moreover, 

ln[l + A]// £—= dp^l (5.13) 

Jo e -e t +p 2 +p 

as A — ► 00. That is why, for simplicity, we take ln[l + A] to evaluate F(A) 
and obtain an approximative radiative correction 

327T 

R C = ae — ln[l + A] (5.14) 
15 

for the binding energy. (Indeed, for A ~ 1 these two functions perfectly 
coincide) 

5.3. Calculating concrete values of a and A. We search for values of 
a and A, that guarantee the error of the self energy Err being smaller than 
the radiative correction Rq- This leads to the condition: 

■ ^e ln[l + A] avr , 
a < Mini -If- 01 L_ ,— V. 5.15 

> 28A 2 [1+A 2 ] 7 a 2 4A 

37T 2 (1 — a) 
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Set A = 1/4 which corresponds to a photon energy mc 2 /2, that is an enor- 
mously high value compared to the binding energy eo- 

Therefore, the condition on the coupling parameter, such that the radia- 
tive correction Rq dominates the error £rr, is given by 

a < 0.85(/?Z) 2 . (5.16) 

Asking for the nuclear charge numbers Z that guarantees enhanced binding 
in the physical case a = 1/137, leads to Z > 13. 

However, in the case of smaller values of A ~ eo, which seems physically 
reasonable, we are in a perfect shape and the error of the self-energy is by 
far negligible compared to the energy shift. 



Appendix A. Auxiliary Lemmas 
LEMMA 1. Let^ eH. Then 

DD\p 2 + HfY 1 D*D*^j < {0\DD[p 2 + Hf}' 1 D* D*\0)\\^\\ 2 

+ const. A^,H f ^) + ||p*|| 2 ). (A.l) 
Proof. We fix an arbitrary photon number n. Recall, 



n+2n+l 

T[D*D*ip n ] n+2 = £ £ G(k 3 ) • G(h) 

v / (n + 2)(n + 1) 



n+2n+l 
3= ' 

X ^(^Ja, . . . . . . . . . ,k n+2 ), (A.2) 



where indicates that the j—th variable is omitted. By permutational 
symmetry we can distinguish between three different terms, 

(Vn, DD[p 2 + H f ]~ 1 D*D*^ = I n + II n + III n . (A.3) 

First, the diagonal part, 

T V- f ■ ^(^Ml, h,-.., fcn+2)| 2 „ 

\,v=l,2 J r l^i=l K 't\ +Li=l \ K i\ 

(A.4) 

If we set Q= \l-J2"=3 +\h + k 2 \ 2 + J2^ 2 \k t \ and b = 2 [l - ££+ 2 h] ■ 
[ki + k 2 ] and use the expansion 

' LM t W (a.5) 



Q-6 Q Q Q Q Q-b Q 
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then we again see that the second term vanishes when integrating over ki,k,2- 

I |2 

Therefore, with Q > \k\ + &2 + \ki\ + \k2\ we arrive at 



I n < 



£ [ 

A,i/=1,2 



dk\dk2 



+ 4 



[G^fcQ-G^fca)] 2 
|fci + fc 2 | 2 + |A;i | + |/c 2 | 



[g^fcQ-G^fca^Q^I + lfcal] 2 

[|fcl + fc2| 2 + |fci| + |fc2|] 2 (|fci| + |fc2|) 



n+2 



|^~y^ *i| 1^(^*3, • • • , k n+2 )\ 2 dldki . . . dk n+2 



i=3 



< (Q^lp 2 + H f ]~ 1 D*D* 



+ const. A||f0„ || 2 . (A.6) 



For convenience we define the operator which can be regarded as 

the norm of D(x), 

\D\(x)=^2 j \G x (k)\e ik - x a x (k)dk (A.7) 



A=l,2 



\D\* denotes the operator adjoint. Similarly, we can define Obvi- 
ously, |GLLl Lemma A. 4] still holds for the "norm" of D and E, namely 



\D\*\D\<-H f , \E\*\E\<^-H f , 

7T 07T 



(A.8) 



that can be proved analogously to (|3.51|) . For the second term, with p 2 + 
Hf > Hf, we evaluate 



Hn<n J2 / 



\G x {k 1 )-G v (k2)\\G x (k l )-G v (k n+ 2)\ 



< 



x \ip n (l,k 3 , . . . , k n+ 2)\\ip n (l, k 2 , ■ ■ ■ ,k n+ i)\dldki . ..dk n+2 
\G(k^ 2 



\k 



^-dfei(|^|,|Z>r|Z>||^ n |) < const. A 2 (A.9) 



Finally, the term where on one side the indices of the created photons 
differ completely from the indices on the other side, 



\G x (h) ■ G»(k 2 )\\G x (k n+1 ) • G"(fc ra+2 )| 

ESN 



III n <n 2 

X,v=l,2 

x \ip n (l, h,---, K+2)\\i>(}, h,---, k n )\dldki . . . dk n+2 
< [W n \,\D\*HJ l/2 \D\*\D\HJ 1/2 \D\W n \) <const.A(^ n ,^ n ), (A.10) 
where we used 



n+2 


71+1 


1/2 


n+2 


1/2 




£n> 

i=l 


i=l 




i=2 




(A.H) 
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the fact that we can write 

'1/2 



[H f ] L,2 Mhh,---,kn) 



t=l 



1/2 



i) n {l,k u ...,k n ), (A.12) 



and the first equation of (|A.8|) . □ 
LEMMA 2. Let f £W and fix an arbitrary photon number n. Then 

y 3 / 2 (o-E*4, n+1 ,\p 2 + H f ]- 1 D*D*^ n ) < 



o 



< const. 



a\A + A 3 )||Vg| 2 + a(A + A 2 )(V>„+i, #/Vwi) 



Proof. Obviously, 



where 



4<(n + l) £ | 



^>n+l, b 2 + H f ]- 1 D*D*?p n ) =I n + II n 
' H\k l )\\G x {k l )\\G v {k 2 )\ 



(A.13) 
(A.14) 



I A; 



■ x 



A,^=l,2 

x \ip n (l,k 3 , . . . , fc n+2 )||^„. + i(Z, fe, • • • , k n+2 )\dldh . . . dk n+2 
<const.A 2 ||Vn|| 1/2 / |G(fe 2 )| ^/NP^S)* 2 

< const^a^A 3 !!^!! 2 + a- 1 / 2 A 2 (^„+i,i//Vn+i)) (A.15) 



For the second term of (|A.14|) we obtain 

\H x (kx)\ \G x (k n+1 )\\G»(k n+2 )\ 



II n <(n + lf/ 2 Yl I 

\ , 1 o " 



It- 1 

X,u=l,2 J 2-<i=l l^l 

x \il> n (lifa, ■ ■ ■ ,k n )\\ip n+1 (l,k 2 , ■ ■ ■ , k n+2 )\dldki . ..dk n+2 

Lff A (fci)||^ n (i, fa,..., fc, 



<(n + l) 3/2 £ / 



A,i/=1,2' 



14, . . . , K n ) 



|G A (fc„+i)||G , '(fc„ + 2)||l/'n+l(i, k 2 , . .. , k n+2 )\ 



1/2 



dldk\ . . . dk n+2 



< \\\E\HJ 1/2 \^ n \\\\\\D\HJ 1/2 \D\\^ n+1 \\\, (A.16) 



which implies the statement of the lemma by use of (|A.8j) . 



□ 
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